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EVALUATION AND EXPANSION OF AN ANALYTICAL MODEL 
FOR FATIGUE OF NOTCHED COMPOSITE LAMINATES 


R, L. Ramkumax, S. V. Kulkarni* 
and R. B. Pipes** 

Materials Sciences Corporation 
SUI4MARY 

This report describes the analytical and experimental 

study performed to expand the existing static and fatigue 
failure analysis . The analytical effort extended the analysis 
to include interlaminar effects , while the experimental effort 
developed methods to obtain basic experimental data required 
as input to the analysis . 

The static failure analysis for notched laminates was 
modified to include interlaminar effects near the notch . The 
interlaminar effects were included by discretizing the laminate 
into laminae in the region of the projection of the notch width 
along the axial direction . Thi s discretization produced addi- 
tional degrees of freedom which were used to analyze axial 
in-plane shear damage or crack in each lamina and axial inter- 
laminar damage or delamination at the laminae interfaces . 

Because a large number of damage zones is possible when elastic- 
perf ectly plastic behavior in in-plane and interlaminar shear 
is considered , the analysis was simplified by assuming a 
secant modulus representation of the nonlinear shear stress- 
shear strain curve. The grovi/th of damage by both the above- 
mentioned approaches is calculated for various realistic 

T-300 Gr/Ep laminates ([0„/±45] , [0/±45/90] ) and compared 

z s s 

with earlier experimental observations. 

Three-dimensional elastic and two-dimensional elastic- 
plastic finite element analyses were also performed for some 
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notched laminates for the purposes of (i) developing an under- 
standing of the behavior of interlaminar stresses around the 
no tch, and ( ii) determining the nature of the 2-D elastic- 
plastic stress distribution adjacent to the notch when the 
maximum attainable stresses or combinations thereof for the 
laminate are governed by a yield criterion. The results 
obtained from task (i) provided a guideline for qualitatively 
assessing the approximate stress analysis procedure. Task 

(ii) determined how appropriate some of the postulates (such 
as the average stress concentration region being a function of 
of laminate anisotropy) of the basic failure model are. 

The experimental effort was directed at several differ- 
ent tasks: (i) determining static/fatigue (S-N) curve data 

for interlaminar shear and normal stresses; (ii) assessing 
the effect of (a) the interaction between in-plane shear and 
transverse compression on the wearout rates for strength and 
stiffness in a [±45]^ laminate subjected to uniaxial loading, 
and (b) the transverse restraint (as provided by 90° layers 

in a [0./90.] laminate) on the failure mode and wearout rates 
1 D s 

for strength and stiffness for a unidirectional layer when 
subjected to compression/compression fatigue; and finally 

(iii) designing the experiments and specimens to obtain the 
required data. 
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INTRODUCTION 


The use of composite materials for primary aerospace 
structural applications requires that they have adequate safety 
margins and lifetimes under all anticipated environments. 

Accordingly , a sound understanding of the fracture and fatigue 
behavior of these materials is needed , and the materials must 
also be well characterized by soundly based experimental data . 
To characterize materials by tests alone is expensive ; analy- 
ses which correlate successfully among data points can greatly 
reduce the number of tests needed and guide the testing 
methods to provide more useful experimental results. 

The development of such analyses was initiated in refer- 
ence 1 in which a methodology was developed for predicting 
the growth of cracks and ultimate failure under fatigue loading 
of fiber composite laminates containing through- the- thickness 
notches or holes . This framework was based on knowledge of 
the static and fatigue behavior of a unidirectional layer of 
the material . Specifically , the analysis assumes that lamina 
fatigue behavior within the notched laminate is predictable 
from fatigue data obtained from tests of unnotched laminae . 

To verify this analysis , and develop a semi -empirical 
capability for the fatigue analysis of notched , f iber-composite 
laminates , the predictions derived from reference 1 were com- 
pared in reference 2 with experimental fatigue data for notched 
[ 02/145 ] ^ Boron/epoxy laminates containing a 0.635 cm diameter 
circular hole . The specific phenomena studied were : 

(1) initiation of fatigue damage , and grov/th as a func- 
tion of load cycles ; 

(2) fatigue life and mode of failure ; and 

(3) residual strength and mode of failure after a pre- 
deterrained number of cycles . 

The residual strength after fatigue loading, and both 
axial and transverse damage growth, were raeasured quantita- 
tively for correlation with the analytical predictions . 



Correlation between analysis and experiment was hampered , 
however, by the following factors: 

(1) de lamination of the 0° surface layers in the region 
of the longitudinal projection of the notch diameter; 

(2) lack of a statisitically significant data base for 
lamina fatigue properties in general ; and , par- 
ticularly , 

(3) absence of lamina axial compression fatigue data . 

The principal deficiency in the analysis appeared to be 
the lack of a capability to predict the growth of the delamina- 
tions . A primary objective of the present study is to correct 
this deficiency . The approach used is to modify the existing 
static failure analysis (ref . 3 ) to permit the effects of 
interlaminar shear stress to be taken into account. This is 
achieved by discretizing the laminate into laminae in the 
vicinity of the notch. Significant features of the analytical 
developments are described in the section "Description of the 
Failure Model , " and the mathematical details are presented in 
Appendix A. The analysis has been incorporated into a com- 
puter program for calculation of damage growth and residual 
strength. 

To assess the degree of realism associated with the 
assumptions in the analysis , the calculated stresses have 
been compared to those obtained from 2-D and 3-D finite ele- 
ment analyses . These comparisons are discussed in the sec- 
tion "Finite Element Analysis of Notched Composite Laminates . " 

An experimental program was conducted to demonstrate 
the feasibility of obtaining the data required as input to 
the analysis . Inadequacies of the existing data base for 
lamina fatigue properties (already noted during the analysis- 
experiment correlation study in reference 2 ) are compounded 
by the need to incorporate consideration of interlaminar 
effects into the fatigue analysis . Exploratory testing to 
develop methods for obtaining lamina properties from simple 
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laminates instead of unidirectional samples therefore has been 
initiated as follows; 

( i ) [ ± 4 5 ] ^ laminate uniaxial compression/ compression 

fatigue tests have been made to determine the nature 
of the interaction between in-plane shear and trans- 
verse compression in a lamina and to compare the 

wearout rate with that of a [±45] uniaxial tension/ 

s 

tension fatigue test; 

( ii ) [0p/90j]g laminate uniaxial compression/compression 

fatigue tests have been made to ascertain the effect 
of transverse restraint (provided by the 90-degree 
laminae ) on the failure mode and the wearout rate 
for a unidirectional lamina in compression/ 
compression fatigue ; 

( iii ) [+ 252 / 90 ] g laminate specimens were used to obtain 

the S-N curve for interlaminar normal tensile stress ; 
and 

(iv) Partially notched [±45 . / 0 . ] laminates were used 

1 [] s 

for determining the behavior of interlaminar shear 
stress in fatigue . 

Highlights of the tests and their implications in the 
light of the analysis as extended herein are discussed in the 
section entitled "Experimental Program . " The test specimen 
des ign is outl ined in Appendix B. 
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FAILURE MODEL FOR NOTCHED COMPOSITE LAMINATES 


INCLUDING INTERLAMINAR EFFECTS 

DISCUSSION OF PREVIOUS MODEL 

The prior analysis (ref . 2 , see fig . 1 ) was based on two 
concepts : (1) Failure in the presence of a notch can result 

from cracks propagating along planes of weakness parallel to 
the fibers or perpendicular to them; (2) The notched laminate 
may be divided into regions having simple , approximate stress 
states . The f irst of these concepts is utilized unchanged in 
this study . The second is extended to take into account 
through- the- thickness effects . 

In reference 2 , the laminate was considered as an assem- 
blage of discrete regions ; namely , the strip containing the 
notch; the overstress region adjacent to this strip; the inter- 
mediate , shear-over stress region; and the undisturbed or aver- 
age stress region . In the thickness direction , the laminate 
was treated as a homogeneous material . Consequently , the pre- 
dictions of interlaminar stresses and the resulting delamina- 
tions were beyond the scope of the analysis . To take into 
account the interlaminar stresses , the laminate was considered 
herein divided into laminae in the local regions where the 
effect of interlaminar stresses is most pronounced . Resultant 
changes in the analysis and its failure-prediction capability 
are discussed in the following section. 

DESCRIPTION OF THE NEW FAILURE MODEL 

Static Failure Analysis 

The static failure model , at the laminate level , is 
divided , as before (ref . 1), into three regions : ( i ) a central 

core region which is the projection of the notch in the loading 
direction; ( ii ) an overstressed region of average stress con- 
centration , adjacent to the core region ; and ( iii ) an average 
stress region (fig. 1 ) . As mentioned in the earlier sections , 
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this model adequately predicts the gross heterogeneous 
behavior of the laminate in the overstressed and average stress 
regions. The stress state in the vicinity of the notchf though, 
indicates the presence of high interlaminar stresses. In some 
laminates, these interlaminar stresses are large enough to 
cause a delamination or "peeling off" of one ply from another 
(ref . 2). This necessitates discretization of the core region 
to the lamina level so that interlaminar stresses may be quan- 
tified and the corresponding failure mode predictive capability 
built into the model (fig, 2). 

The predictable failure modes in the analysis comprise : 

( i ) an axial in-plane crack in any layer in the core region; 

( ii ) an interlaminar debond between any two layers in the core 
region; and (iii) a transverse crack across the laminate, 
through the specimen thickness , normal to the loading direc- 
tion. The failure modes are thus predicted at the lamina 
level in contrast to the gross laminate failure modes of the 
earlier model (fig. 3) . 

The governing equations for the laminates that exhibit 
the above failure modes form five coupled equations (Appendix 
A) , The core region is assumed to have three different plies 
above the midplane . The mathematical complexity introduced 
by the actual nonlinear shear behavior of the individual 
laminae is overcome by assuming an elast ic -per feet ly plastic 
behavior and an ultimate strain criterion for crack initiation . 
Alternatively, an elastic , secant modulus behavior may also 
be assumed . The latter approach reduces the computational 
efforts by a considerable amount. The words "defect" and 
"damage" are used to refer to in-plane and interlaminar plastic 
regions , and in-plane crack and interlaminar delaminations , in 
general. This facilitates a general discussion on both the 
elastic-plastic approach and the elastic , secant modulus 
approach . The governing equations and their solutions for 
regions in the laminate v/ith different damage zones are given 
in Appendix A. 
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The elastic solution for the notched laminate (Appendix A) 
estimates the maximum in-plane and interlaminar shear stresses 
at the notch . A comparison of these stresses , normalized with 
respect to the corresponding yield values , determines the type 
and number of defects that initiate and grow simultaneously in 
the initial damage zone and the corresponding magnitude of the 
load. Two or more defects that occur in succession due to a 
small change in the applied loading are assumed to initiate and 
grow together to enhance the predictive capability of the model . 
The growth of the initial damage zone is predicted through the 
solutions to the modified governing equations (Appendix A) . 

These solutions are compatible with the elastic solutions at 
the boundary of the two regions (fig. 4 ) . The normalized maxi- 
mum shear stresses and shear strains at the notch are computed 
for each increment of the first damage zone size to predict 
the initiation of the second damage zone and the corresponding 
load . 

The second damage zone in the notched laminate may have 
additional inelastic layers and interfaces , or may have a crack 
or debond in a layer or interface that was inelastic in the 
first damage zone . The growth of the second damage zone size 
w-ith the applied loading involves solutions for three regions 
that are continuous at two boundaries (fig. 4 ) . This procedure 
can be extended to many successive damage zones till the lam- 
inate fails . A transverse failure is marked by a crack across 
the specimen , through the thickness , normal to the loading 
direction . An axial failure is determined by an unstable 
growth of a crack/debond with the applied loading . 

Figure 5 shows the various possible sequences of defects 
in a general notched laminate. The branching of defects is a 
strong function of the laminate material and geometric proper- 
ties . The solutions for any damage zone depend on the solu- 
tions for the regions that preceded it . Consequently , general- 
ization of the solution procedure is very tedious when more 
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than two damage zones exist in the laminate. However , this 
is alleviated by eliminating a few damage sequences that are 
unlikely to occur . It is assumed that , before too many damage 
zones appear , the load will be large enough to cause a trans- 
verse failure or one of the earlier defects will precipitate 
an axial failure . This is where the simpl ic ity in the applica- 
tion of the elastic , secant modulus approach is obvious . The 
elimination of the inelastic regions leaves only cracks and 
del ami nations to be accounted for . 

In many laminate constructions, a degenerate region is 
created whenever the applied loading reaches a critical value. 
Degeneracy sets in when one or more layers in the core region 
are surrounded by stress-free surfaces and/or inelastic sur- 
faces . In the elastic-plastic approach, one or more degenerate 
layers have a linear variation of axial strain in the degenerate 
region. The degenerate layers in the elastic, secant modulus 
approach are stress free and are essentially "carried" by 
adjacent load-carrying regions. In a degenerate region, the 
solutions for the degenerate layers are uncoupled from the 
five basic governing equations. The solutions change for dif- 
ferent types of degeneracies (Appendix A) . Hence, if all types 
of laminate behavior are to be included in the solution pro- 
cedure, programming of the solution becomes tedious. 

Fatigue Failure Analysis 

A "mechanistic wearout" concept used in earlier studies 
(ref . 2 ) underlines the basic philosophy of the laminate 
fatigue behavior . On fatigue loading, material property degra- 
dation is predominant in the vicinity of the notch due to 
stress concentration effects . An experimental characteriza- 
tion of the lamina fatigue data helps estimate the degrada- 
tion in the laminate properties due to cyclic loading . The 
material property degradation on fatigue loading , when incor- 
porated into the static failure model described earlier, could 
lead to fatigue failure modes and strengths that are very 
different from the static predictions . 
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Figure 6 shows a flow chart of the program developed to 
predict the failure modes and strengths of notched laminates 
sub j ected to cyclic loading . 


DISCUSSION OF RESULTS 


The approximate failure analysis for notched laminates 
which considers the effects of interlaminar stresses was 
utilized to investigate the damage growth in some T-300/5208 
laminates . Table 1 lists the various data sets utilized in 
the failure analysis . The trial data are used to illustrate 
the extent of the predictive capability of the analysis when 
laminate failure in the transverse direction is prevented . 

Tables 2 and 3 indicate the results of the analysis for 
T-300/5208 [O^Z-^B/O] ^ and [±45/02/0] ^ laminates, respectively 
for different values of interlaminar boundary layer zone size 
"d" and the trial data. In the table, the distance "d" is 
the interlaminar boundary layer region which assumes two 
values s (i) half the lamina thickness (typical for Boron/ 
epoxy); and (ii) one-tenth the lamina thickness (typical for 
Graph! te/epoxy) . Figure 7 shows a typical variation of the 
in-plane displacement as obtained in reference 4. The figure 
explains the meaning of the parameter "d". In the present 
analysis, the axial displacements in the various layers are 
piecewise uniform, and hence, discontinuous at the interfaces. 
The ratio of the displacement discontinuity to "d" - the trans 
verse shear strain - should match the slope of the continuous 
displacement curve at the interface obtained from a higher 
order analysis . For the [ 02/±45/0] ^ laminate , the in-plane 
plastic zone is initiated almost simultaneous ly in the O 2 and 
0 laminae ^^ 0 ) • Also , the next damage zone - a plastic 

zone - appears in the ±45 laminae (a^y2^®^ because the inter- 
laminar shear stress (t ,„) is not large enough to form an 
interlaminar inelastic region. Thus , the approximation that 
the axial , in-plane plastic zone grows uniformly through-the- 
thickness (ref . 3 ) may reflect the actual behavior for this 
laminate . After both the plastic zones have grown over a 
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certain distance, an in-plane crack is initiated in the 0 ^ 

and 0 layers . Note, however , that the interlaminar shear 
stresses are high and may eventually give rise to a delamina- 
tion . In tables 2 and 3 , the results for the [0„/±45/0] 

2 s 

laminate are essentially the same as for the [±45/0„/0] 

z s 

laminate except that the magnitude of the interlaminar shear 

stress (T is increased by about 25%. 

xzl2 

An important conclusion which can be drawn from the above 
results is that the number of combinations of different damage 
zones (inelastic regions as well as cracks) is large (fig. 5) . 
The factors which affect this number are the material proper- 
ties, stacking sequence, and notch size. An automated numer- 
ical scheme to "track" the various possible dam.age growth 
combinations is a tedious task . Also , the solution procedure 
expands with the initiation of plastic zones . Thus the analy- 
sis is further complicated . Consequently , the analytical 
procedure was modified to eliminate the need to consider 
plastic zones by considering linear approximations of the 
nonlinear in-plane and inter lairiinar shear stress-strain curves 
as shown in figure 8 . Such an analysis would , however , still 
consider the growth of in-plane shear cracks and interlaminar 
debonds , but the number of combinations of the different 
damage zones is drastically reduced . The various damage growth 
combinations are illustrated in figure 5 . 

Figures 9 and 10 illustrate the progressive damage growth 
for [0„/±45/0] and [0/±45/90] T-300/5208 Graphite/epoxy 

Z 3 S 

laminates, respectively, based on the trial data in table 1. 

The force represents the applied laminate force and "a" 
indicates the cumulative sum of the damage zones. In figure 
9, two inelastic in-plane axial damage zones are predicted 
before an axial crack in the 0° layers appears. Figure 10 
illustrates a similar pattern of damage growth, except that 
no in-plane crack occurs in the 90° layer and an interlaminar 
inelastic zone is initiated at a higher value of the applied 
force . Figure 11 shows the effect of notch size on the damage 
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growth in a [±45/0/0] laminate for the trial data input. The 

s 

[±45/0/0] notation is used for a [±45/0„] laminate, which 
s ^ s 

is treated by the analysis as a "±45° layer" and two 0° layers. 
If the ±45 layup is treated as two layers (a +45 layer and 
a -45 layer) in the analysis, the Poisson effect causing an 
increase in the shear modulus will not be accounted for. 

Figures 9 to 11 show the capability of the approximate analysis 
procedure to predict damage regions on a lamina level . How- 
ever, as has been noted earlier, the analysis is not capable 
of predicting all combinations of damage growth because of a 
significant increase in the complexity of the computation 
procedure . 

The variation of failure load with notch size is studied 
in figure 12 for a [±45/0/0]^ laminate for the actual data 
input. The following conclusions can be drawn from the figure: 

(i) Damage zone size increases with increase in notch 
size . 

(ii) In-plane damage occurs only in the 0° layers. 

(iii) Laminate fails by transverse crack propagation. 

(iv) At failure, the maximum in-plane shear stress in 
the ±45 laminae is higher for larger notch sizes 
than that for smaller notch sizes, thus indicating 
the increased probability of through-the-thickness 
axial crack propagation for larger notch sizes. 

(v) The transverse failure values predicted by the 

secant modulus approach are close to those predic- 
ted by the elastic-plastic analysis. 

Figure 13 illustrates the results for a [0/0/ ±45]^ 
laminate with the allowable interlaminar shear stress value 
reduced to account for the presence of interlaminar normal 
tension stress at the 0/±45 interface. The obvious result 
is the initiation of an interlaminar inelastic zone followed 
by in-plane plasticity in the 0° layers. 
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In figure 14 the results for a [0/±45/90] laminate for 

s 

different notch sizes are shown for the actual data . Since 
this laminate is more isotropic than the [0 „/ ±45 ] ^ laminate , 
the inelastic damage zones are smaller . Also , failure occurs 
by transverse crack propagation. 

Figure 14 is an illustration of the predictive capability 
of the secant modulus analysis . As shown in the previous 
figure, the [0/±45/90] laminate fails in the transverse 
direction before any crack appears when the secant modulus 
approach was used. In figure 15, the damages are made to 
grow by preventing a transverse failure. Table 4 shows the 
effect of the interlaminar boundary layer thickness on the 
damage growth for the same laminate. 

Figure 16 illustrates an interesting interlaminar damage 
growth in a [0/0/±45]^ laminate for a reduced value of inter- 
laminar stress that causes delamination. The reduction 
accounts for the effect of interlaminar normal tensile stress. 
A high value of the applied stress at which transverse 

failure occurs, is chosen to prevent that failure mode. It 
is seen that the delamination that initiates at the 0/+45 
interface grows in an unstable fashion to failure. The 0° 
plies thus separate away from the ±45 plies in the notch 
region. This result was observed experimentally (fig. 17) 
in the first phase of the program (ref. 2). This qualitative 
agreement between the failure prediction of the analytical 
model and the experimental observation justifies the approach 
taken in building the model. 

In the results discussed above , the ultimate failure 
stress for the ±45 laminate has been obtained from "netting 
analysis . " The netting analysis predictions may be sig- 
nificantly higher than first-ply failure predictions because 
fiber failure within a ply is required in "netting analysis . " 
Also, a through- the- thickness axial crack cannot propagate 
in a ±45 laminate unless the f ibers fail . Laminate failure 
stress corresponding to the fracture of ±45 layer fibers is 
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synonymous to the netting analysis failure stress . Hence , the 
upper-bound predictions of the netting analysis may be realis- 
tic when considering the incipient conditions required for 
the initiation and growth of through- the- thickness axial cracks . 
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FINITE ELEMENT ANALYSIS OF 
NOTCHED COMPOSITE LAMINATES 

The failure predictive capability of the analytical model 
was based on a few simplifying assumptions. Justification of 

these assumptions is essential to extend the model application 
to any general laminate. 

A three-dimensional finite element analysis was carried 
out to predict the in-plane and interlaminar stresses in a 
notched laminate. The primary objectives of the 3-D finite 
element analysis were (i) to substantiate the assumptions 
made in the approximate model and ( ii ) to provide guidance 
for any modification . This elastic analysis aimed at estab- 
lishing a qualitative agreement with the model stress pre- 
dictions , and also establishing the predominance of the axial 
normal and interlaminar shear stresses over the other stresses 
A two-dimensional, el as tic -perfectly plastic finite element 
analysis was carried out to verify the inelastic damage growth 
predictions of the analytical model for various notched 
laminates . 

THREE-DIMENSIONAL ELASTIC ANALYSIS 

The present analytical modeling philosophy considers the 
physics of the problem without accounting for the detailed 
stress distribution . Thus , only those stress fields (and the 
appropriate equations of equilibrium) are considered which 
have a predominant effect on the observed failure modes . The 
driving force for this approximate approach has been the com- 
plexity of the three-dimensional nonlinear stress analysis and 
the associated failure prediction. However , the use of simpli 
fied models (which neglect the effects of some stress fields ) 
can be justified only if there exists qualitative , and to a 
limited extent , quantitative correlation with the more exact 
stress analyses . To that end , a three-dimensional finite 
elem.ent anlays is of notched laminates was carried out to check 
the assumptions made in the approximate stress analysis pro- 
cedure . 



The computer code SAP- IV was used to carry out the finite 
element analysis (ref . 5) . An octant of the notched laminate 
was considered and each layer above the reference midplane 
was discretized into the mesh form shown in figure 18 . A 
minimum of eight and a maximum of twenty-one nodes were chosen 
to describe the general three-dimensional isoparametric ele- 
ment used in the analysis. The element was used to represent 
orthotropic , elastic media. Laminate analysis was used to 
apply the appropriate loads in each layer so that a far-field 
uniform strain state existed in the laminate. 

In figure 19, obtained from the approximate analysis, 

at the 0/±45 interface, is compared with the variation 

from the 3-D finite element analysis (the average value at 
the face centroids) . The comparison is made at two axial 
locations. As shown in figure 19, the average stresses from 
the approximate analysis compare qualitatively with the 3-D 
finite element results. 

Figure 20 shows the axial normal stress variation in the 

load-carrying direction agreeing qualitatively with the analyt 

cal prediction. Figure 21 shows the variation of the same 

stress, a , in the elements just above the x axis. A com- 

y Y 

parison with the applied stress in each layer (fig. 18) gives 
a measure of the stress concentration in each of the three 
layers of the [0/+45/-45]^ laminate. These stress concentra- 
tion values show good agreement with the value predicted by 
the analytical model. 

The number of integration points to be used in the numeri 
cal evaluation of integrals over volumes is referred to as the 
integration order in SAP-IV. Different orders of integration 
may be chosen in the plane of the element and through its 
thickness . The orders of integration chosen in the finite 
element analysis affect the cost of the run and the magnitude 
of the results . Figures 22 and 23 show that for two different 
orders of integration ( 2 and 4 ) , the variations with y at 

the 0/45 interface for a [0/±45/-45 ] ^ laminate have good quali 
tative agreement along with a quantitative difference . 
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The approximate analysis , discussed in the previous chap- 
ter, treats notches of any geometry as equivalent slit notches 
since the notched strength for composite laminates is not 
strongly dependent on the notch shape (ref . 3 ) . The 3-D 
finite element analysis was used to verify this . Figure 24 
compares the shear stress variations for a circular hole and 
a slit notch of the same width at the ±45/-45 interface of 
the [0/±45/-45]g laminate . The results show a slight differ- 
ence only in the vicinity of the notch , justifying the notch 
geometry assumption in the simplified analytical model . 

The 3-D finite element results indicate that the variation 
of in-plane stresses (a and o ) justifies the assumption 

L yy Xy 

of an average stress concentration region and a localized 

region of shear stress transfer. The finite element solutions 

also justify the approximate analysis assumption of a uniform 

value of (independent of the x-coordinate but a function 

of the y-coordinate) in the notched region and a zero value 

in the other regions of the laminate. The non-zero value of 

a predicted by the 3-D finite element analysis in the region 
yz 

away from the notch is of little significance because it does 
not contribute to the precipitation of failure. 

2-D ELASTIC-PLASTIC FINITE ELEMENT ANALYSIS 

As an additional task conducted during the present study, 
a 2-D elastic-plastic finite element analysis was performed 
for some notched laminates . Because of the consideration of 
el as tic- perfectly plastic behavior in shear in the basic 
fracture model , the elastic-plastic finite element analysis 
was used to investigate the damage growth and stress distribu- 
tion near the notch and to ascertain the effect of laminate 
anisotropy on damage propagation . The PLANE finite element 
program (developed at Grumman , ref . 6 ) was used for this 
purpose . 

Figure 25 illustrates the finite element mesh and the 
loading for a laminate containing a slit notch . Input values 
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of moduli and strength (ref . 7 ) for various laminates are shown 
in table 5 . Note that the failure stresses satisfy the stabil- 
ity condition for the Hill ' s yield criterion (ref . 8) . Results 
of the finite element analysis are presented for a [0/±45] ^ 
laminate in figures 26 and 27. The notch blunting effect due 
to plasticity is evident in both the figures. The figures also 
show the change in the notch blunting effect on varying the 
applied load. Subsequently, the plastic region growth is inves- 
tigated for laminates with various anisotropies. The critical 
value of the loads (a multiple for the applied load vectors) at 
which plasticity is initiated for the [0], [0/90]^ and [02/±45]^ 

T-300/5208 laminates is tabulated in table 6. As expected, 
the load is lower for the [0] and [0/90] ^ laminates. In 
figures 28 through 30 are shown the elements which have yielded 
for all the three laminates for similar load levels. It is 
evident that the yield zones are most prominent for the [0/90] ^ 
and [0] laminates, in that order. However, it has been ob- 
served in experiments with Boron/epoxy laminates that the 
axial inelastic damage zone is more pronounced in the [0] 

laminates than in the [0/90] laminate (ref. 9) . This is 

s 

also predicted by the static failure model, when the effec- 
tive in-plane shear modulus is modified for laminates contain- 
ing 90° layers in order to reflect the local bending effects 
of 90° fibers in the vicinity of the notch (ref. 3) . The 2-D 
finite element analysis does not account for these changes. 
Hence, it appears that the elastic-plastic finite element 
analysis is not capable of predicting accurately the growth 
of damage zone for the [0] laminate. Table 7 shows the varia- 
tion of an approximate measure of stress concentration factor 
with increasing load for different laminates. The stress con- 
centration shows a decrease with increasing load for the 

[0/90] and [0„/±45] laminate, while, for the [0] laminate, 
s ^ s 

it increases . This increase is contrary to the well-known 
notch blunting effect when there is a growth of the inelastic 
region in the axial direction . 
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EXPERIMENTAL PROGRAM 


The experimental tasks conducted during this prograin have 
been defined by the observations made during the previous 
analysis-experiment correlation study and the expansion of the 
basic failure model to include interlaminar effects. The 
tests are of an exploratory nature and were performed for the 
specific purpose of defining the requirements of basic data 
to be utilized as input in the analysis. 

The thrust of the earlier experimental program (ref. 2) 
was to provide Boron/epoxy lamina static/fatigue data to be 
used as input into the fatigue analysis and to generate static/ 
fatigue data for notched Boron/epoxy laminates to be utilized 
in the fatigue failure model. However, in the previous programs 
several items were not investigated . 

Consequently , there is a need for uniaxial compression 
fatigue data , even for a tens ion/ tens ion fatigue loading of a 
notched laminate . To that end, [0] coupon compression tests 
were performed . However , since failure in a [0] longitudinal 
compression test is accelerated by the presence of longitudinal 
cracks , the residual strength/stiffness data obtained may not 
be representative of a unidirectional layer in a multi-layered 
laminate , because of the absence of the restraining effect 
provided by the adjacent layers . Hence , utilization of the 
unidirectional compression data alone to predict laminate 
fatigue response is questionable . Alterations in failure modes 
for static compression have been observed in reference 10 on 
providing a nominal lateral reinforcement . For these reasons , 
it may be imperative to obtain the compressive fatigue data 
in the presence of transverse restraint. Indeed , just as the 
[ ±4 5 ] g specimens are used for shear data , another simple 
laminate may be required for other layer properties . For axial 


compression , 


[ 0 ./90 . ] 
r 3 s 


laminate provided 


Comparison of [0./90.] and [ 0 ] 

1 3 ^ 

test results will indicate the 


laminate 
effect of 


compression fatigue 
lateral constraint on 


residual properties . (Note that the residual properties for the 
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0° layers in the [0./90.] specimen have to be extracted by 

1 D 

filtering out the contribution of 90° layers) . 

Combined stress effects upon wearout have not been evalu- 
ated experimentally and have been modeled approximately in 
the fatigue model. Hence, initial assessment of the magnitude 
of this problem was obtained from compression fatigue tests 
on [±45]g laminates. These results can be compared with those 
obtained from tension fatigue tests on similar laminates. Thus 
axial in-plane shear combined with transverse tensile stress 
will be compared with the same shear and transverse compressive 
stress . 

Important parameters in the analysis which considers debond 
propagation are the interlaminar strengths (shear and normal 
tension) . An estimate of the interlaminar shear strength may 
be obtained from the short beam test. Alternatively, the 
interlaminar shear strength was obtained from a [±45/0] ^ lam- 
inate with a slit (transverse to the fibers) in the 0° layers. 
Also, as suggested in reference 11, the interlaminar (a^) 
normal tension strength was determined by utilizing a [±252/90] ^ 

laminate. The [±25/90] laminate is optimized to yield a 

^ s 

maximum value of a^. 

The various tests and the specimen geometries are summar- 
ized in table 8. 

MATERIALS FABRICATION AND QUALITY CONTROL 

Fabrication of the composite laminates was accomplished 
with contemporary vacuum bag autoclave techniques with the 
proper cure cycle as supplied by the material producer. In 
order to evaluate process conditions for adequacy, lamina 
characterization tests were performed to determine fiber ten- 
sion and in-plane shear properties. These properties are 
summarized in table 9. Review of the lamina characterization 
tests indicated that the fabricated T-300/5208 composite 
possessed an average fiber tensile strength of 1668 MPa, a 
Young's modulus of 141 GPa , and an in-plane shear modulus of 
5 GPa. These properties are indicative of high-quality 
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laminate fabrication (ref . 7) and suggest that the processing 
conditions were adequate. 

Each composite panel was subjected to ultrasonic inspec- 
tion after fabrication. In this way , both defects could be 
isolated from test specimens. A typical "C" scan of a 32-ply 
[±45]^ laminate is shown in figure 31. The "C" scan reveals 
anomolies near the , "raw" panel edges associated with variations 
in resin content or thickness due to edge effects. These 
regions of the panel were removed prior to specimen fabrica- 
tion. Interior regions are revealed as uniform in both thick- 
ness and resin content. 

All panels were fabricated from a single batch of the 
T-300/5208 prepreg. The panels were fabricated during the 
initial phase of the program and stored in dessicant for 
future use. Specimen fabrication was accomplished by first 
cutting the composite panel into sublaminates. The sublaminate 
was then mounted in a fixture for bonding tabs. After tabs 
were bonded to the sublaminates, test specimens were sawed 
to the proper geometry with a precision diamond saw facility. 
Each specimen was then engraved for identification as to 
laminate, sublaminate, and position within the sublaminate. 

Next, specimens were instrumented with electrical resistance 
strain gages and stored for subsequent testing. 

TEST METHODS 

The test methods developed in this study include methods 
for determination of 0° compression (constrained and uncon- 
strained) fatigue properties, the in-plane shear fatigue proper- 
ties in the presence of a compressive transverse stress, inter- 
laminar shear fatigue properties , and interlaminar normal 
fatigue properties of the T-300/5208 material systems. 

Compression Tests 

The test method chosen for evaluation of the 0° compres- 
sion and in-plane shear fatigue properties was a modified 
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IITRI wedge action grip test fixture shown in figure 32 . The 
IITRI compression fixture consists of wedge action friction 
grips mounted in support blocks , with connecting rods and 
linear bearings which yield a concentric load path . While 
the IITRI fixture was originally designed for a narrow test 
specimen of 6.35 mm in width with an unsupported length of 
6.35 - 12.7 inches, the fixture was modified to accept speci- 
men widths up to 19.05 mm. In addition, specimen lateral 
deflection restraint was achieved by special roller supports 
attached to the wedge action grips as shown in figure 33. 

In the evaluation of 0° compression strength of the 
T-300/5208 composite system, the test specimen geometry 
strongly influences test results. In addition to unsupported 
specimen length, the specimen width was found to affect com- 
pressive strength. For example, the average strength for a 
specimen geometry of 12.7 mm width and 28.575 mm unsupported 
length was 957 MPa, while the average strength for the 6.35 
mm width specimen was 1134 MPa. By reducing the unsupported 
length of the 6.35 mm width specimen to 15.875 mm, the average 
compressive strength achieved was 1347 MPa with a corresponding 
modulus of 133 GPa . These results compared favorably with the 
tensile properties (1668 MPa and 141 GPa) and, therefore, the 
last specimen geometry was chosen for fatigue testing. The 
test specimen's unsupported length is important in the deter- 
mination of compressive strength due to restraint of axial 
cracking rather than column instability about the axis of mini- 
mum moment of inertia. Indeed, the failure mode of the 0° 
compression test specimens was axial splitting followed by 
micro- instability of the segments formed by the axial cracking . 
It would appear, therefore, that the compressive properties 
developed in this effort are uniquely a function of the test 
specimen geometry. The actual compressive strength of the 
composite specimen is bonded by the Euler critical load for the 
total specimen and the strength achieved when axial cracking 
is fully constrained . 
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To investigate the influence of axial cracking upon the 
compressive strength of the composite, the [0/90/0], 
laminate was found to exceed a "rule of mixtures" prediction 
of strength (935 MPa) and thereby illustrate the synergistic 

influence of the 90° constraint layers . In order to minimize 
the influence of edge effects upon the behavior of the 
[ 0 / 90 / 0 ] 2g test specimen , a specimen width of 19.05 mm was 
chosen . Since it was shown earlier that specimen width 
influenced the compressive strength of the [ 0 ] test speci- 
men, it should be expected that the large width would lead 
to conservative estimates of the [0/90/0 ] specimen strength . 
Hence , the difference between observed strength and "rule of 
mixtures " estimate ( 7 % ) could have been greater had the speci- 
men width been smaller and the associated deleterious influ- 
ence of edge effects been absent. In order to evaluate the 
influence of unsupported length upon compressive strength of 
the [ 0 / 90 / 0 ] 2g laminate, specimen tests with unsupported 
lengths of 12.7 mm and 9.525 mm were conducted. The differ- 
ence in average strengths for the two test lengths was less 
than one percent. It would appear, therefore, that the 
lateral support system was effective over such a variation 
in unsupported length and acted to insure compressive failure . 
It should be noted , however , that the typical stress- strain 
response revealed some bad eccentricity as exhibited by dif- 
ferences in apparent moduli between stress-strain responses 
of gages on opposite sides of the spec imen . Failure commonly 
corresponded to an abrupt increase in strain on one s ide of 
the specimen and a corresponding decrease on the other side . 
This behavior would tend to imply an instability induced 
failure . However , the lateral constraint at the center of 
the specimen should preclude gross instability. Hence, the 
failure of the [ 0/90/0 ] 2^ specimen is considered to be ini- 
tiated by a local instability of the outer 0 ° layers . 

In order to investigate the in-plane shear fatigue proper 
ties of the T-300/5208 material system in the presence of 
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a compressive transverse normal stress, [±45,,] laminate 
specimens were loaded in compression. Since the compression 
response of the [±45]^ laminate is of nonlinear softening 
character, it was necessary to employ a 32-ply laminate in 
order to preclude instability induced failure. A specimen 
width of 19.05 mm was chosen in order to minimize the influ- 
ence of edge effects upon laminate response. By introducing 
a teflon layer between specimen and tab adjacent to the test 
section, the test section could be effectively increased while 
maintaining the lateral constraint of the tab. In order to 
assess the influence of the transverse normal stress upon 
static in-plane shear strength, it is informative to compare 
the tensile and compressive strengths of the [±45]^ laminate. 
The average static compressive strength observed was 201 MPa 
with an average modulus of 18.9 GPa . The average static 
tensile strength for the identical specimen geometry was 170 
MPa with a modulus of 18.8 GPa. The latter data show good 
agreement with the tensile response of the 8-ply [±45]^ 
laminate of 25.4 mm width and 152.4 mm test section length. 
Average tensile strength and modulus for this specimen geometry 
were found to be 150 MPa and 19.2 GPa, respectively. Hence, 
i-t would appear that there is a strong dependence of shear 
strength upon the sign of the transverse normal stress. 

The laminate specimen was subjected to compres- 

sion fatigue employing the modified IITRI test fixture shown 
in figure 33. Both fatigue life and residual strength data 
were developed with this test method. Residual strength test 
specimens were instrumented with longitudinal and transverse 
strain gages for determination of residual modulus and 
Poisson's ratio. 


Tension Tests 

Tensile test methods were developed and evaluated for 
determination of the interlaminar shear fatigue, and inter- 
laminar normal fatigue properties. The interlaminar shear 
fatigue test specimens consisted of a [±45/0^]^ laminate 


24 



wherein the 0° layers were discontinuous at the center of the 
specimen. Load carried by the 0° layers was transferred 
through interlaminar shear . Hence , failure was initiated at 
the discontinuity of the 0° layers in an interlaminar shear 
mode. Specimens were subjected to a fatigue loading (R = 0.1, 
30 Hz) and the initiation of cracking detected visually. 

Static samples were instrumented with acoustic emission sen- 
sors for detection of interlaminar failure. Specimens were 
not tested beyond the initiation of damage. 

The interlaminar normal stress fatigue properties were 
evaluated with the [±252/90]^ tensile coupon. Due to the 
large difference in Poisson's ratios between the 90° and ±25° 
layers, failure is initiated at the specimen edge in an inter- 
laminar tensile mode. Fatigue specimens were subjected to 
tensile fatigue loadings of R = 0.1 at a rate of 30 Hz. Both 
acoustic emission and visual detection techniques were employed 
to establish failure initiation. 

DISCUSSION OF RESULTS 

As noted above, the objective of the uniaxial compression 

static/fatigue tests of [±45]^ Graphite/epoxy laminates is to 

determine the effect of interaction between in-plane shear and 

transverse stress (transverse tension plus in-plane shear for 

a tension test of [±45] laminate vs. transverse compression 

s 

and in-plane shear for the compression test) . To that end, 

static tension and compression tests were conducted for the 

[±45] laminate and the results are tabulated in table 10. 
s 

It is seen that the average compression strength is higher than 

the average tension strength, although the moduli are almost 

identical. This is an expected result and is confirmed by 

analytical predictions. Thus the [±45]^ laminate cannot be 

used realiably for predicting shear strength because the mode 

of transverse stress has a noticeable effect on it. 

Table 11 and figure 34 show the summary of the T-300/5208 

[±45^,1 laminate compression fatigue data . Note that there 
J_ o s 
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is only a slight wearout of strength and stiffness. This is 
contrary to the significant wearout observed for the tension 
fatigue tests . Note also that the scatter is rather high 
for the lifetime data at S = 0.67, Figures 35 and 36 
illustrate the failure modes in static tension, and static 
compression after fatigue, respectively. The failure in ten- 
sion occurs primarily due to ply separation which is precipi- 
tated by transverse tension failure in the ±45 laminae. Micro- 
buckling at the lamina level is the predominant failure mode 
for compression. An interesting failure mode is observed in 
figure 37 for a [±45]^ laminate subjected to compression. The 
surface ply has buckled in figure 37 resulting in further loss 
of stiffness and subsequent failure. 

The static/fatigue data for interlaminar normal stress 

a as obtained from a [±25„/90] laminate are shown in tables 
z 2. s 

12 and 13 and in figure 38. In table 12, the initiation of 
failure at the free edges was monitored by acoustic emission 
and from the load-displacement chart. In figure 38, the 
scatter is rather large for S = 0.6. However, this is to be 
expected. The slope of the S-N curve does indicate a wearout 
rate that is consistent with a matrix-dominated behavior. The 
edge failure mode for this specimen is illustrated in figure 39. 

The static/fatigue test data for the [± 45 / 02 ] ^ partially 
notched laminates are given in tables 14 and 15 and shown in 
figure 40. Failure is initiated in the interlaminar shear 
mode and this is evident from the C-scan in figure 41 and the 
photomicrograph in figure 42. Referring to the derivations 
in Appendix B and to table Bl , the analytical prediction based 
on an interlaminar shear strength of 68.95 MPa is substantiated 
by this data. Interestingly enough, the scatter for S = 0.8 
is rather significant while for S=0.5, the data points are 
rather close. Because of the notch blunting effect during 
fatigue, the residual strength for the specimen at S = 0.6 
is higher than the average static strength. Figures 43 and 
44 illustrate the static and residual strength, and fatigue 



failures, respectively . The fatigue failures , hov/ever , are 
visible only under magnification. 

The static compression results for the [0] and [0/90/0] 
laminates are given in tables 16 and 17 . The compression 
strength for the [0] laminate is consistent with the values 
reported in the literature. The modulus is also comparable 
to the [0] tension modulus as reported earlier in table 9. 

The compression strength and modulus for the [0/90/0] 2 ^ 
laminate is comparable to the values obtained by multiplying 
the corresponding quantities for the [0] laminate by 0.66 
(percentage of 0 degree laminae in the laminate) . 

Lifetime and residual strength data for the [0] laminate 
are tabulated in table 18 and plotted in figure 45. An impor- 
tant feature of the [0] fatigue tests is that for most of the 
specimens, failures occurred by the separation of tabs because 
of the high stress levels. Hence only two fatigue failures 
were obtained. When a specimen experienced separation of tabs, 
it was removed from the fixture, the tabs were rebonded, and 
the specimen was tested statically to failure. Hence the 
result is a residual strength data point. Because of these 
difficulties, an S-N curve for the [0] laminate could not be 
obtained. Note also that the residual strength actually shows 
an increase in strength. This appears to be paradoxical. A 
probable explanation for this inconsistency is the scatter in 
the static strength data. The compression fatigue failures 
for the [0] laminate are illustrated in figure 46. The failure 
sequence is splitting in the axial direction, instability, 
and brooming. 

The [0/90/0] 2g compression fatigue test results are shown 

in table 19 and figures 47 and 48. Once again, the data points 

for S = 0.8 and 0.67 sometimes differ by two or three decades. 

While the residual strength for the [0] shows an increase for 

S = 0.5, that for the [0/90/0] 2 ^ laminates registers a small 

decrease. Note, however, that the failure modes for [0] and 

[0/90/0] „ laminates are different. In figure 49, the residual 
z s 


27 



strength test failure suggests the buckling of the surface 0 
degree layer on an elastic foundation followed by overall 
instability and subsequent fracture . The fatigue failures in 
f igure 49 for S = 0.8 occurred generally inside the tabs , 
suggesting that some of the data points in figure 47 are pre- 
mature failures. 
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CONCLUDING REJ'4ARKS 


The basic failure model for predicting through-the- 
thickness crack growth and ultimate failure of notched com- 
posite laminates (ref. 3) has been expanded to include the 
effects of interlaminar stresses in the vicinity of the notch. 
The modification was accomplished by locally discretizing the 
laminate into laminae in the region of the projection of the 
notch diameter along the direction of the loading. In-plane 
and interlaminar damage (inelastic region or crack) growth 
was calculated for various laminates ([0^/±45]^, [02/±45/0]^ 

and [0/±45/90]^) and stacking sequences. 

Because of the possibility of a large combination of 
in-plane and interlaminar damage zones and the numerical com- 
plexity associated with "tracking" them, the capability of 
the failure model is restricted to the prediction of only a 
limited number of combinations of these regions. In order 
to simplify the analysis, a secant modulus approximation to 
the nonlinear in-plane and interlaminar shear stress-strain 
curves was considered in lieu of the elastic, perfectly 
plastic approximation. This eliminated the need to consider 
inelastic regions and damage was predicted as a crack in a 
layer or at the interface. This enhanced the capability of 
the analysis to predict axial crack growth in an unstable 
fashion or until failure. 

Three-dimensional finite element analysis was performed 

in order to develop an understanding of the nature of the 

stress distribution around a hole in [0/±45] and [±45/0] 

s s 

laminates. Ths SAP-IV computer program, was utilized for 
this purpose. The results obtained from the finite element 
analysis were compared with that of the approximate stress 
analysis procedure with the objective of qualitatively assess- 
ing the kinematic assumptions. 

As an additional task, a 2-D elastic-plastic finite 
element analysis was conducted for some notched laminates. 


29 



Because of the consideration of inelastic effects in shear in 
the basic failure model, it appeared desirable to conduct such 
an analysis for the purpose of investigating the nature of the 
damage growth and the stress distribution in the vicinity of 
the notch, and to ascertain the effect of laminate anisotropy 
on damage propagation. 

The experimental program basically consisted of the fol- 
lowing exploratory tests: 

(i) [±45]^ laminate uniaxial compression/compression 

fatigue tests to determine the nature of the interac- 
tion between in-plane shear and transverse compres- 
sion in a lamina and to compare the wearout rates 
with that of a [±45]^ uniaxial tension/tension 
fatigue test. The results indicated that there is 
no noticeable wearout, thus suggesting that the 
sign of the transverse stress has an important effect. 

(ii) [ 0 ./ 90 .] laminate uniaxial compression/compression 
1 j s 

fatigue tests to ascertain the effect of transverse 
restraint (provided by the 90 degree laminae) on the 
failure mode and the wearout rate for a unidirec- 
tional lamina in compression/compression fatigue. 

The data scatter for the [0/90/0]^ laminate 
fatigue tests was rather large for S = 0.8 and 
specimen failures were experienced outside the gage 
section. Problems in the [0] laminate tests were 
premature failures along the bond between the speci- 
men and the tabs, and an increase in the residual 
strength after fatigue. The former could. be attrib- 
uted to high stress level while the latter could be 
a result of a scatter in the static strength data. 
Because of these inconsistencies, an appropriate 
evaluation of the [0/90/0] 2 g laminate to obtain 
unidirectional compression fatigue data could not 
be made . 
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(iii) [±252/90]^ laminate specinien to obtain the S-N 
curve for interlaminar normal tensile stress. 

(iv) Partially notched [± 45 ^/ 0 ^]^ laminate for deter- 
mining the behavior of interlaminar shear stress 
in fatigue. 

The current analysis includes the effects of damage growth 
due to both in-plane and interlaminar shear stress behavior. 

The transverse stresses in the notched laminate are ignored 
in the analysis in comparison to the axial stresses. This 
assumption was validated by finite element analyses and experi- 
mental observations . The experimental observations on a 
notched [ 02/±45 ] ^ laminate were predicted qualitatively by 
this analytical model using an elastic , secant modulus approach . 

The shortcoming of the failure prediction methodology is 
the extensive effort needed to automate the numerical proce- 
dure to account for all possible damage sequences in a general 
laminate . A recommended alternative , for future research , 
would be the utilization of these defined damage and failure 
modes in a finite element analysis . This could provide a 
better understanding of the stress state around the notch and 
have the additional advantage of incorporating all the stresses 
for general in-plane loading into the failure prediction 
methodology . 
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APPENDIX A 

FAILURE MODEL FOR NOTCHED COMPOSITE LAMINATES 
WITH INTERLAMINAR EFFECTS 


Referring to figures 1 and 2 in the main text, the govern- 
ing equilibrium equations for the various regions in the notched 
laminate may be written as shown belov/. 


ELASTIC REGION 

The basic approach utilized in formulating the present model 
is to discretize the laminate into laminae only in the region of 
the projection of the notch width. The remaining portion of the 
homogeneous laminate is divided into an average stress concentra- 
tion region and a uniform stress region as in reference 2 . It is 
also assumed that the symmetric laminate has a maximum of six 
laminae. Equilibrium of each ply in the central core region and 
the adjacent overstressed region yields the following equations: 


d^u-, 2G . 

01 _ xyl ( 


dx 


E ,aa,, ^^01 ^1^ 

xl 0 


*^xzl2 , \ _ n 

E .h^d ^ 01 02^ ° 

xl 1 1 


(Al) 


d^u 


dx 


02 

2“ 


2G 


G 


xy2 , ■> , ^xzl2 , 

E .aa^ ^^02 " ^1^ E h d ^"oi ^02' 

x2 0 x2 2 1 


(u - u ) = 0 
E „h„d„ 02 03 

x2 2 2 


(A2; 


d^u„_ 2G ^ 
0 3 xy3 


2G 


(u„^ - u, ) + 


xz23 


dx 


E ^aa„ ' 03 1' E„^h^d 


x3 0 


x3 3 2 


^^02 ^ 03 ^ ° 
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d u. 


dx"’ 


2G T h ^ 


E 


xul 0 


2 , ^^01 

a„h 


, ^'^xv 2^2 . 

"'"2 ■ ^^^02 
^xul^O^ 


Ui) 


+ 


G ^h_ 

xy3 3 


E 


xul^O 


2 , ^^^03 

a„h 


Ui) - 


xyl 2 

E ,a^ 
xul 0 




“ 2 > 


(A4) 


Note that the inplane shear stress transfer region is of 
width "aQ" for all layers and the interlaminar shear stress trans- 
fer is assumed to take place over a thickness "d” . Axial force 
equilibrium in the laminate requires that: 


du 


2ahn E T 
1 xl 


01 


dx 


+ 2ah„E ^ 
2 x 2 


dx 3^x3 dx 


+ 2a hE - 5 — 

0 xul dx 


du-, dUp 

+ 2a„hE „ = F 

2 xu 2 dx 


(A5) 


Integrating equation (A5) with respect to x yields; 


2«2'>^xu2“2<*> = Fx + - 2ah^Ej^j_u„px) 


^^■’2‘=x2"02<^' - ®>’3®X3“03'=‘> ' 


where is a constant of integration. 

Substituting equation (A 6 ) into (A4), the governing equations 
can be written in a matrix form as: 


d^{u) 


dx" 


[B]^fuTg - (Fx + Cj^q) {d} 


where {u}^ J^u^^ Uq2 


T 


(A7) 

(A 8 ) 
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and the subscript or a superscipt "e" represents the elastic 

region . 

The general solution of equation (A7 ) is : 

-m-, X m, X -m X m x 

{u} = (C,,e + C_e ) {U} + (C e ^ + C e ^ 

e 11 16 1/ t / 


-m^x m^x 

+ (C^3e + C^ge ){U} + (C^^e 


-m,x 


m.x ,,, 

+ C^ge ^ ){U} + (Fx + C^q){S}^ 


(A9) 


where {S} 


[B]g ^{Dl 


(AlO) 


m 


^,m 2 ,m 3 and are the eigenvalues of the matrix [Bj^, 


and are the corresponding 

eigenvectors 

and *“10 ' *“11 ' *“12 ' constants. 

Finiteness of the stresses for large x is introduced by 
setting ^ y _ 1 ' *“18 *“19 zero. Hence 


lu-iX { ~[ \ m^x /o\ ^ ^ ' 

{u} = ^ {u} + C,„e lu) + C, .e {u> ' 

e 11 12 


'13 


+ Cj^^e {U} + (Fx + C^q) {S}^ 


(All) 


Five boundary conditions are required to solve for the five con- 
stants. If the applied load is such that the stresses in the 
laminate are in the elastic range , these conditions are the zero 
stress condition on the notch surface (x = 0) and symmetry 


34 



conditions in the and regions at the same location . That 
X s ^ 





and u® (0) = u® (0) = 0 


(A12) 


INELASTIC REGION (S) (0<x < a) 

Three inplane and two interlaminar inelastic shear stress 
transfer regions are postulated in the discretized region. If, 
for example, the top layer (layer number 1) has an inplane 
inelastic region of length a (see fig. 1 ) , equations (Al) and 
(A4) are modified in accordance with the elastic-perf ectly plastic 
shear behavior assumption (fig. 6 ): 


d u 


01 


dx' 


xzl2 

E , h- d^ 
xl 1 1 


(u 


01 


u 


02 ' 


"xl- 


(A13) 


.2 

d Ui 


dx' 


+ 


2G „h„ 
xy2 2 

E ^ a^h 
xul 0 


(u 


02 


G 0^0 

- U,) . -2Zl_| 


E ah 
xul 0 


u 


03 


Ui) 




xul 0 


2 ^^1 


u. 


2 T h-, 
xyl 1 

E . a^h 
xul 0 


(A14) 


Equations (A2) , (A3) and (A5) remain unaltered. 

For an inelastic region at an interface (see fig. 1 ), a 
similar modification is introduced . Assuming an interlaminar 
inelastic region between layers 1 and 2, for example , equations 
(Al) and (A2 ) are modified as shown below: 
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( AX 5 ) 


U- ^ T 2 VJ 

01 


dx' 


(u 

Exiaao 01 


. ^xzl2 

“i> “ e;7f- 


ilj>2 _ ^^xy2 , 

,2 E '^02 

dx x2 0 


G 


u^) 


xz23 


u. 




^03^ 


'xzl2 

E „h^ 
x2 2 


(A16) 


Equations (A3), (A4) and (A5) remain the same. 

The governing equations for a region of the laminate where 
inelasticity is present may then be written as: 


d ( u j- 


P 


dx 


[ B ] { u } 

P P 


(Fx + C^) {D} + {R} 


(A17 ) 


where {R} represents the inelastic contribution to the nonhomo- 
geneous terra. The solution for this set of equations is: 


-rn X nvx , , ^ -m,x m x , ^ ^ 

(u) = (C,e ^ -t- C,e ^ ){Ul^ ^ T (C„e ° + C.,e ^ ){u}^ 

p 1 6 2 / 

-m X m x ,3. -m X 

+ (C,,e ' + C„e ' + (C.e 

3 o 4 


rn X 

+ Cge ® + (Fx + C 5 ){S}p + {T} 


(Ais: 


where {T} = -[B]^^{R] and {S} = -[B] ^{Dl (Al9) 

P P P 

and the subscript or the superscript "p" represents the inelastic 

region . 

The eigenvalues m.(i = 5,8) and the corresponding eigenvectors 
/ i ) ^ 

{U}' of [B]., for the inelastic region are different from those 

for the elastic region due to the changes in the [B] matrix. 
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Coordinate x for each region is measured from the origin of that 
region . 

Consider now an inelastic region in the laminate , extending 
over a length of a from the notch. Equations (All ) and (Al8 ) 
are the displacement solutions for the elastic and inelastic 
regions. The boundary conditions for the problem are; 



In addition to the above 14 boundary conditions, the inelas- 
ticity condition for the appropriate layer/interface at the junc- 
tion of the two regions must be imposed for a unique F - a 
relationship . 

REGION WITH AN INPLANE OR AN INTERLAMINAR CRACK (0 < x < C) 

The governing equations for a region with a crack are 
obtained by imposing the stress-free conditions on the cracked 
inplane or interlaminar surface . For an inplane crack in the 
top layer , for example , equations (Al) and (A4 ) are modified to ; 


37 



(A21) 


,2 

d u 


dx‘ 


01 

9 


xzl2 

E ^h^d. 

xl 1 J 


iU 


01 


02 ' 


2 

d Ut 2G 

I + xy2 2 , 

,2 2, ^ 02 
dx E , a^h 
xul 0 


, , '^xy3^3 , 

“l> + ; 2T <‘>03 


Ui) 


^xyl2 


E , 
xul 0 


(u- 


u- 


Equations (A2), (A3) and (A5) remain unchanged. 

Iff on the other hand, there is a debond between layers 1 
and 2, equations (Al) and (A2) take the form: 


^ ^01 
j 2 

dx 


2G 


xyl 


E^lUaQ 


(u 


01 


u^) 


0 


d u 


02 


dx 


2G . 
xy2 

E „aaf, 
x2 0 


(u 


02 




xz2 3 

E „hod„ 
x2 2 2 


, u 


02 




Equations (A3), (A4) and (A5) remain unchanged. 

For either type of crack, the governing equations take the 
form of equation (A7) , and the solutions, the form of equation 
(A9) . Introducing new constants, 


d^{u] 


dx" 


[B] {u) = (Fx + C' I-) (Dl 

c 25 


(A22) 


(A23) 


(A24) 


(A25) 
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“m„x 




' 27 ' 


m _ 

+ C,-j ) 

z b 

{U} 

(1) 

' ' -I- ( n 

^ 22 

l0X)^U}^2) 

+ 

-m, , X 

/ J — L 

^^23® 


“^IqX 


’^11^ (3) 

+ C2ge ){U} 


-m X m X 

+ (C 24 S + {U}^ ‘ + (Fx + C^^){S}^ (A26) 


where {S} = -[B] ^ D (A27) 

c c 

and the superscript or subscript "c" represents the cracked 
region . 

The boundary conditions for a problem with a cracked region, 
an inelastic region, and an elastic region could be written in a 
manner similar to equations (A20) . For a unique F - a - C 
relationship, the criterion for crack initiation will also have 
to be imposed at the crack tip. 

DEGENERATE CASES 

The various regions in the laminate as described above are 
determined by the sequence in which defects are likely to be 
initiated with an increase in the applied load. However, more 
than one defect could initiate and grow simultaneously. For 
some widely used laminate constructions, certain sequences of 
defect formation could cause the nature of the above solutions 
to change. Such a situation will arise (i) when two or more 
inelastic zones (inplane/interlaminar) appear simultaneously, or 
(ii) when one or more laminae act as non-load carrying members 
undergoing rigid body displacements. The first type of degeneracy 
occurs for elastic-plastic behavior in inplane/interlaminar shear, 
and the degenerate layer/layers have a linear variation of axial 
strain. The second type of degeneracy is more likely to occur 
in an elastic analysis with a shear behavior governed by a secant 
modulus approach. The degenerate layer/layers in this case are 
essentially "carried" by the load-carrying adjacent regions. 
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Simultaneous Inelastic Inplane and Interlaminar Shear Stress 

Transfer for a Lamina 

Consider, for example, the first type of degeneracy when 
inelastic zones have initiated simultaneously in the surface 
layer 1 (inplane shear mode) and at the interface between layers 
1 and 2 (interlaminar shear mode). Equation (Al) takes the form: 


d u 


01 


dx' 


2 t , 
xyl 

E , a 
xl 


+ = R 


(A28) 


Integrating this equation with respect to x twice, one obtains 


Uqi(x) 


R 


1 2 


+ C3QX + 


(A29) 


If the boundary of this region for x = 0 is the notch surface, 
then C 3 Q vanishes. Equation (A2 9) then becomes: 


^01 


(x) 




'31 


(A30) 


Equation (A2) is then modified to: 


d^u 


02 


2G 


xy2 


dx 


E „aa. 
x2 0 


(u 


02 


Uf) - 


xz23 

E ^h„d 
x 2 2 2 


(u 


02 


^03^ 


T 


E 


xzl 2 

x2^2 


and equation (A4) assumes the form: 


(A31) 
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d u, 2G „h„ 
i , xyz z 
^ + 

, a^h 
xul 0 


(u 


dx E 


02 


G ^ ^ o 

Ui) + (Uq 3 

^xul^O^ 


Ui) 


'xyl2 


E ia„ 
xul 0 


2 <"l 


“2> = 


2t hi, 
xyl 1 

‘=xul®o‘' 


(A32) 


Upon utilizing the solution in equation (A30) , equation (A6) 
becomes 


2t t 

2a„hE „u„(x) = Fx + - ah^E x^ + 


2 xu2 2 


'35 


1x1 E „ a E , 

xl xl 1 


2®‘’2'^x 2"02‘^> - ^‘'3'=x3“03 


2ao‘'Exul“l<>‘) 


or 


2a_hE „u„ (x) 
2 xu2 2 


Fx + C r- ~ ( 2 T I h 4- T ,„a)x 

35 xyl 1 xzl2 


2ah„E ^u.„ (x) - ah^E (x) 

2 x2 02 3 x3 03 


2a„hE u^ (x) 
0 xul 1 


(A33) 


Substituting equation (A3 3 ) into equation (A32) , and considering 
equations (A31) and (A3 ) , the governing equations for the 
degenerate region are written as: 


d^{ u} 


d 


dx 


[BlJul, = (Fx + C^^){d}^ +{Dx2},x + f R} 

da 35 d d 


d 


(A34) 


where { u} 


d 


.^02 "03 


i T 


(A35 : 


41 



and the subscript or superscript "d" refers to the degenerate 
solution ( s ) . The solution to equation (A34 ) is : 


-m, m X ,,, -m, .x 

(u)^ = (C 22® ){U} + 


' 36 ^ 


' 33 ' 


+ C3^e ){u}^ ’ + (€340 ^ + C3ge ^ 


where {F} = -[B] , {Dx2}, 

d d d 


+ {F}.x + (Fx + + {T}. 

d 35 d d 


-1 


(A36) 

(A37) 


{S}^ = -[B] {D}^ 


(A38) 


{T}^ = -[B]"^({R} - 2{F}) . 


(A39) 


The boundary conditions for a laminate with a degenerate 
region are written in the same form as before. 

Simultaneous Inplane and Interlaminar Cracking 

for a Lamina 

If, instead of the elastic-plastic approach, an elastic 
secant modulus approach is chosen to predict the "peeling off" 
of the surface lamina (e) in a laminate, the above solutions could 
be modified accordingly. Considering the same example presented 
earlier, the surface ply degeneracy of the second type changes 
equations (A30) through (A33) by eliminating all the terms with 
yield stresses . The governing equations then reduce to the 
form: 


d {u} 


d 


dx 




(A40) 
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The solution to equation (A4 0 ) is : 


-m, m X , , . .X 

{u}^ = ^^ 3^0 + C3ge ){U} ^ + (0336 


m, .X -m, pX m x , ^ . 

+ €3^6 ){U} + (*-34® "^^ 38 ^ )iU] 


+ (Fx + C 33 ) {S}^ 


(A4i: 


It is seen that a single ply degeneracy reduces the number 
of coupled equations by one, as reflected by the matrix 
If a group of layers has a degenerate solution in a region of 
the laminate, the matrix reduces in size further. The equa- 

tions and the solutions for these cases could then be written 
as above. 


Simultaneous Inplane Plasticity in All Layers 

An extreme case of inplane inelasticity in all the layers 
over a region (0 < x < a) of the laminate is not uncommon. In 
that region, the core region (see fig. 1 ) experiences a plastic 

flow. Mathematically, this uncouples the core displacement equa- 
tions from those for the u^^ and u^ displacements. Equation (a 4) 
for this region is: 


ah. 


dx 


’xyl 2 


E .a,, 
xul 0 


2 ^^1 


u. 


* ^^xy2^2 ^ ’^xy3*’3* ^ 




or 


d u 


dx 


®xul^0 


b. 


""xul^O^ 


(A42; 
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(A43) 


where + '^xyS^S- 

Equation (A 6 ) becomes: 

2a„hE „u„ = Fx + C.„ - b^ x^ - 2a hE u (A44) 

2 xu2 2 40 1 0 xui i 

The last two equations can be solved simultaneously for u^ (x) 
and ^ 2 (x) in the degenerate region of the laminate. 

The core displacements in this region are governed by the 
plastic flow equations. Each layer has an inplane shear force 
that is proportional to a, and an interlaminar shear strain that 
maintains the value at x = a. Therefore, the total axial force 
in each layer, or the axial strain in that layer, varies 
linearly with respect to x. That is, 

du . 

— 3 -^ = E.x , i = 1,3 (A45) 

dx 1 


where E -ah.,E- 
xl 11 


2 T h, 
xyl 1 




(A46: 


E „ah„E„ 
x 2 2 2 


^^xy2^2 




(A47; 


and E _ah_E_ = 2x h - 2x (a)a (A48) 

x3 33 xy3 3 xz23 

The core displacement solutions for the degenerate region 
are dependent on the displacement solutions of the adjacent region 
(along the x-axis) through the ’’■xz23^“^ terms. 

In other words. 
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Tx^ 12 (“) /degenerate region = 

T ^ ^ ( a) /degenerate region = 
xz23 


Tx 2 i 2 ^ /adjacent region 

T ^ ^ ( 0 ) /ad j acent region 
xz23 


(A49) 


When an elastic analysis with a secant modulus shear behavior 
is carried out, a region (0 < x < C) of through-the-thickness 
crack might be predicted. For this region, equations (A42) and 
(A44) become: 


d^u. 


dx' 


xyl2 

E a^ 
xul 0 


(Ui - U2) 


(A50) 


2a„hE „u„ 
2 xu2 2 


Fx + C 


40 


- 2a hE ^u^ 
0 xul 1 


(A51) 


respectively. The above equations yield solutions for u^ (x) and 
U 2 (x) for the region 0 < x < ?. 

The core, in this case, moves as a rigid body, and the dis- 
placements are those of the adjacent region along the x-axis . 

UQi/degenerate region = Uq^ ( 0)/adjacent region, i = 1,3 (A52) 
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APPENDIX B 
TEST SPECIMEN DESIGN 


The objective of the experimental program is primarily to 
perform exploratory tests for determining basic data required as 
input into the failure model. After considering the deficiencies 
in the data during the previous analytical/experimental correla- 
tion study (ref. 2 ) and the needs for the modified basic failure 

model with interlaminar effects, the following tests were identi- 
fied : 


(i) [±45]^ laminate uniaxial compression/compression fatigue 

tests to determine the nature of the interaction between inplane 

shear and transverse compression in a lamina and to compare the 

wearout rate with that of a [±45] uniaxial tension/tension 

s 

fatigue test ; 

(ii) [0./90.] laminate unaxial compression/compression 
1 3 s 

fatigue tests to ascertain the effect of transverse restraint 
(provided by the 90 degree laminae) on the failure mode and the 
wearout rate for a unidirectional lamina in compression/compres- 
sion fatigue; 

(iii) [±25„/90] laminate specimen to obtain the S-N curve 
z s 

for interlaminar normal tensile stress and 

(iv) Partially notched [±45^/0 laminate for determining 
the behavior of interlaminar shear stress in fatigue. 

With regard to the design of all specimens subjected to com- 
pression, dimensions were chosen such that failure would not 
occur prematurely due to instability. For the [0./90-] speci- 
mens, the stacking sequence was [ 0 / 90 / 0 ] This would avoid 
any interlaminar edge failures due to compression loading. The 
design of the [±25^/90]^ laminate was taken from reference 
This particular stacking sequence results in a high value of the 
interlaminar normal tensile stress when the specimen is subjected 
to tension loading. This ensures edge failure due to interlam- 
inar normal stress. Finally, a partially notched [±45./0.] 

1 3 s 
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laminate , with the notch in the 0 degree layers , was designed 

to investigate the interlaminar shear behavior , 

Figure Bl illustrates a partially notched [±45 , /O . ] lamin- 

1 J s 

ate with the notch in the 0 degree layers across the width of 
the specimen. Consideration of equilibrium in the x-direction 
results in the following differential equations; 


^ dx^ 


u. 


•u 


+ 


1 G 


xz = 0 


(Bl) 


,2 

dx^ 


’^2 '^1 G 


xz = 0 


(B2) 


where ' d ' is the width of the interlaminar shear transfer zone. 
Subscripts 1 and 2 correspond to the 0 degree and ±45 degree 
plies , respectively. Also, force equilibrium at any cross sec- 
tion of the laminate results in the following expression for 
the displacement U 2 : 


u. 


rFX 

^2w 






(B3) 


E^h u (0) 

where the constant = - 2 ^ applied 

laminate force. Upon substitution of equation (B3) in equa- 
tion (Bl) and the subsequent solution of the resulting dif- 
ferential equation, one obtains 

-» fv X 

u^ = C 2 ® + Ax + B ^ ■£ 
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where is a constant of integration. 


2G 

X2 r , X 1 


(B5 


! 


FG 


A 


xz 


WE^E^h^h^a d 


and B 


■2C, G 
1 xz 


E^E2h^h2da/ 


(B6) 


The boundary condition ensuring the finiteness of the force for 
large x has already been considered by eliminating exponential 
terms of the type Consideration of the boundary condition 

du 

^ = 0 at X = 0 yields 

C 2 = A/a (B7) 

and consequently, 

-FG 

C - 2 (B8) 

^ 2Wa(E„h_da -G ) 

2 2 X z 


An expression for the displacement u^ can be obtained from 
equations (B4), (B6) and (B3) . Having known the displacements 

and , the interlaminar shear stress can be calculated as 


T 

xz 


(X) 


u (x) -u (x) 

G [-i— — 
xz a 


(B9) 


The objective of the above analysis is to determine the 
percentage of the 0 degree layers in the [±45^/0^]^ laminate 
such that failure occurs due to the interlaminar stress x 

X z 

in the vicinity of the notch and not due to overstress in 
the ±45° layers . Also failure should not occur in the 
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unnotched portion of the laminate . 

Table (Bl) shows the results of the analysis. Assuming 
that the ultimate value of the interlaminar stress x xs 

^ -X. Z/ 

68.94 MN/M^ (10 ksi) , the 40%/60% [±45/0] ^ laminate appears 

to be the logical choice . A problem associated with the stack- 
ing sequence for this laminate is that for tension loading, the 
interlaminar normal stress at the free edge would be tensile . 

The effect of a tensile a can be alleviated by choosing a 

z 

wider specimen (say, 3.81 cm or 1-1/2 in) . 
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Table 1. 



G (0°) in GPa 
xy 

G (±45°) in GPa 
xy 

G (90°) in GPa 
xy 

G in GPa 
xy 

E (0°) in GPa 

X 

E (±45°) in GPa 

X 

E (90° ) in GPa 

X 

a (mm) 


a (mm) 
o 

W (mm) 


d=d^=d 2 (mm) 


T in MPa 


T ' (0° ) in MPa 

xy 

T (±45° ) in MPa 
xy 


35.935 


6.550 


3.275 


137.896 


22.394 


10.342 


6.350 


2.540 


40.640 


0.070 


68.948 


68.948 


448.162 


I 


35.935 


6.550 


5.516 


137.896 


22.394 


10.342 


6.35 


3.33 


38.10 


0.01 


68.94 


68.94 


723.954 



35,935 


3.283 


5.516 


137.896 


22 .394 


10.342 


6.350 


3.332 


38.100 


0.014 


68 , 948 


68 . 948 


723.954 


xy (90°) in MPa 


68.948 


68.948 


Orn (0/±45/90) in MPa 
j. s 

(0/0/±45) in MPa 
i s 


4481.620 


4481.620 


482.636 


689.480 


Y (0°) 

'xy 


0.021 


0.021 


(± 45 °) 


0.0206 


Y (90°) 

xy 


1.29 


1.6914 


0.025 








Table 2 . 


Stresses in a [0„/±45/0] T-300/52 


For Different Interlaminar Boundar 

d (mm) 




■ xyl (MPa) 


xy2 (MPa) 


xy3 (MPa) 


xzl2 (MPa) 


xz23 (MPa) 


SCF (Stress 

Concentration 

Factor) 


0.070 


68.95 


366.78 

67.51 

19.15 

-5.97 


2.113 


0.035 


68.95 


372 . 21 
68 . 18 
20.22 
- 5.98 


2.116 


a , o>0; a ^>0 
xyl, 3 xy2— 


"'’xyl (MPa) 

68 . 95 

68.95 

^xy2 (MPa) 

448.16 

448.16 

^xy3 (MPa) 

68 . 95 

68.95 

^xzl2 (MPa) 

23.48 

24.44 

^^z23 (MPa) 

-7.32 

-7.23 

“xyl,3("™^ 

0.52 

0.48 

SCF 

2.1098 

2.113 

“xyl, 3^°' “xy2^°' 

C n -,>0 

^xyl, 3- 


"xyl, 3'"“' 

1.93 

1,92 

“xy2 

1.30 

1.36 

i 

"xyl 

68 , 95 

i 

1 68.95 

T (MPa) 

xy2 

448.16 

1 448.16 

T - (MPa) 
xy3 

68 . 95 

68.95 

T , „ (MPa) 

xzl2 

36.13 

38.07 

T (MPa) 

xz23 

-11.25 

-11.27 


0.014 

68.95 

375.57 

68.60 

21.17 

-5.93 

2.117 


68.95 

448.16 

68.95 

25.34 

-7.11 

0.45 

2.114 


1.92 
1.40 

68 . 95 
448.16 

68.95 
40.11 
-11.40 
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Table 3 . 


Stresses in a f±45/0„/0] T-300/520S Laminate For 

2 S 

Different Interlaminar Boundary Layer Thicknesses 
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NOTES FOR TABLE 2 AND TABLE 3 

TABLE 2: 

a^, width of the average stress concentration region 
= 1.27 mm. 

. In-plane shear yield stresses = 68.95 MPa for a 0® layer 

448.16 MPa for a ±45° layer 

. Interlaminar shear yield stresses = 68.95 MPa 

. Subscripts 1, 2 and 3 refer to the O 2 , ±45 and 0 laminae, 
respectively 

Stress Concentration Factor (SCF) is the normal stress 
in the average stress concentration region divided by 
the gross applied stress 

TABLE 3 : 

. Subscripts 1, 2 and 3 refer to the ±45, 0^ and F laminae, 
respectively . 
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Table 4. Damage Growth in a [0/±45/90] T~300/5208 

s 

Laminate for Different Interlaminar Boundary 
Layer Thicknesses (Secant Modulus Approach)* 


d = 0.001'^ 

0 cm . 

d = 0.00279 cm. 

d = 0. 

0 06 9 9 cm . 

Force 

r 

(cm. ) 

Force 

y 

S 

( cm . ) 

F O jT C G 

r 

(cm. ) 

(KN) 

(KN) 

25. 94 

0 . 0 * * 

25.96 

0.0** 

26. 03 

Co 0 * * 

27.25 

. 0254 

27.23 

, 0254 

27.19 

. 0254 

2 7.60 

.0635 

21 . 90 

. 1016 

27 . 87 

. 1016 

2 8. Id 

.12 7 0 

28. 69 

. 2032 

28 , 67 

.2032 

2 9.04 

.2540 

29.41 

. 304 8 

2 9 . 3 9 

. 3048 

2 9.88 

.3810 

30.77 

. 508 0 

3 0.78 

, 5080 


*The applied, laminate stress ^ at which transverse crack 

propates from the notch tip is assumed to be very large. 

**rnplane crack initiates in the 0° and ±45° layers. 











Table 7. 


Variation of Stress Concentration 
Factor for Various Laminates as a 

Function of Applied Load 


Laminate 

Applied Load 

Stress in Element 12 
Applied Load 

[0] 

20.82 

1.87 


29.15 

1.92 


33 . 32 

2.0 

[0/90]^ 

17.54 

1.77 


21.04 

1.77 


28.06 

1.65 

[0„/±45]^ 
2 . S 

30.85 

2.98 


43.19 

2.43 


49 . 37 

2.14 
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■able 8.. ¥arious Tests and Specimen Geometries 



* Length x Width x 'Thickiiess 
* *3 -Static; F-Fatigue 
■‘•'*jnsupported Length ofc Specimen 





Table 9. NARMCO T-300/5208 Graphite Epoxy 
Unidirectional Laminate Tension 
Properties 


Specimen 

E (GPa) 

X 

V 

xy 

0 ^ (MPa) 

X 

t 

c 

X 

L2 

137.9 

.293 

1606 . 7 

11,000 

R3 

137.9 

. 313 

1661.1 

11,700 

LI 

143.4 

.237 

1669.2 

10,930 

R4 

146.9 

. 310 

1735 . 3 

11,220 

Average 

141.6 

.288 

1668.1 

11,210 
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Table 10 . [±45,^] T-300/5208 Laminate Static Test 

16 s 

Data for Uniaxial Tension and Compression 


Specimen 

Type of Test 

a (MPa) 

X 

(GPa) 

V 

xy 

IBRl 

Uniaxial 

Compression 

202 . 8 

~ 

- 

1BR5 

It 

206 . 3 

19 . 8 

. 7645 

1BR2 

n 

192.5 

18.0 

.7565 

ICCl 

n 

203.5 

- 

- 

1CC2 

II 

202.0 

- 

- 

ICC 3 

II 

202.9 


_ 

' 

Average 201.7 

18.9 

.7605 

1BR6 

Uniaxial 

Tension 

171.7 

19.4 

0 .6795 

6AR1 

II 

155.1 

19 . 0 

.8000 

6AR2 

11 

150.3 

19.5 

.7400 

6AR3 

II 

143.8 

18 . 4 

..7500 


Average 155.2 

19 . 0 

.7420 
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Table 11. T-300/5208 Laminate Fatigue 

Compression Data (R = 0.1, f = 30 Hz) 



Residual Compress 

ion Test 

Specimen 

S 

N (cycles ) 

(MPa) 

X 

E^.(GPa) 

V 

xy 

1BR3 

.80 

12,490 




IBLI 

. 80 

4,450 




1BL3 

.80 

7,200 




Avg . 






1AC3 

. 67 

2,146,700 




1BR4 

. 67 

15,450 




1CL3 

. 67 

55,650 




1AC4 

. 67 

1,308,130 




ICLl 

. 67 

500,000 

199.9 

17.58 

. 781 

1BL6 

.67 

500,000 

195.4 

16.55 

.808 

1AL5 

. 67 

500,000 

201.7 

17.44 

. 866 




Avg . 199.0 

17.17 

. 818 

1CL4 

.67 

150,000 

200.4 

— 

— 

1AC5 

. 67 

150,000 

197.7 

17.65 

.788 

1BL4 

.67 

150,000 

192.8 

14 . 55 

.827 




Avg. 197.0 

16.13 

.808 

1CL6 . 

. 50 

10,000,000 

201.5 

18.68 

.819 

1AL6 

. 50 

10,000 , 000 

200.5 

17.65 

. 871 

1AL4 

. 50 

10,000,000 

198 . 7 

17.51 

.777 




Avg . 200.2 

17.93 

. 822 

1AL2 

. 50 

1,000,000 

198.0 

17.10 

. 827 

1AC6 

. 50 

1,000,000 

198.3 

20.20 

.775 

1BL5 

. 50 

1,000,000 

194.2 

18 . 00 

.700 


Avg . 196.8 

1 b . 4 1 

.767 




















Table 13 . 


Tension/Tension Lifetime 

Test Data for [±25/90] 

2 s 

Laminates (f = 30Hz, R = 


Fatigue 

T-300/5208 

0 . 1 ) 


Specimen 

S 

N (cycles) 

2BL2 

0.8 

4510 

2AR2 

CO 

o 

5550 

2AL1 

0 . 8 

312 0 

2AL2 

0 . 7 

7600 

2BR3 

0.7 

8400 

2BL1 

0 . 7 

8000 

2AR3 

0.6 

48,660 

2BR2 

0.6 

205,550 

2BR1 

0 . 6 

794,050 

2CR1 

0.6 

2,061,580 

2DL1 

0.6 

155,600 


64 








Table 14 . 


Static Test Data for 

[± 45 / 0 _] Partially Notched 
3 s 

D-300/5208 Lamii 


1 

! 

Type of Test 

0 ^ (MPa) 

Specimen 

X 

3AR2 

! 

Uniaxial 

Tension 

104.1 

3AL2 

fi 

77.2 

3AL3 

tf 

92.4 



Avg . 91.01 








Table 15 . Tension/Tension Lifetime Fatigue 

Data for [±45/0^] Partially Notched 

3 s 

T-300/5208 Laminate (R = 0.1, f - 30Hz) 


Specimen 

S 

N (cycles) 

0 ^ (MPa) 
X 

3AR1 

0 . 8 

530,450 


3AL1 

0 . 8 

82,120 


3C2 

0 . 8 

8,910 


3BL3 

0 . 7 

504,210 


3BL2 

0 . 7 

2,957,100 


3BR3 

0 . 7 

617,740 


3BR2 

0 . 6 

4,397,520 


3BL1 

0.6 

5x10^ (runout) 

113* 

3C1 

0.6 

5x10^ (runout) 

- 


*Residual 
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Table 17 


Uniaxial Compression Static 
Test Results for [0/90/0] 

s 

T-300/5208 Laminate 


Specimen 

(MPa) 

X 

E (GPa) 

!X 

V 

xy 

5BC1 

1061.5 

87.6 

0.066 

SAC 5 

1038 . 4 

97.9 

0.091 

5AL5 

867.0 

91.0 

0.056 

■SAC 4 

1000 . 8 

93.1 

0.105 

SBLl 

913.2 

91 . 7 

0 .067 

SAC6 

1172 . 1 


— 


1008.7 

92.4 

0.077 
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Table 18 . Uniaxial Compression/Compression 

Lifetime and Residual Strength/Stiffness 
Data for [0] T-300/5208 Laminate 
(f = 30Hz, R = 10) 



Residual 

Residual 
Modulus 
Ex (GPa) 

Specimen 

S 

N, Cycles 

Strength 

(MPa) 

X 

7D14 

CO 

o 

1,420 

1418* 


7D15 

M 

1,470 

FF 


7.D16 

M 

850 

1146* 


7D25 

«l 

4,240 

1217* 


7D26 

II 

4,110 

1656* 


7D2 7 

M 

3,350 

1628* 


7D11 

0.67 

10^ 

1274* 


7D12 

H 

94,510 

1243* 


7D13 

II 

4,677,110 

1192* 


7D22 

II 

10,300 

1772* 


7D23 

II 

10,240 

1487* 


msm 

II 

7,780 

1573* 



0.5 

108,200 

FF 



II 

5 X 10^ 

1288 

- 


II 

5 X 10^ 

1513 

140 


II 

5 X 10^ 

1578 

137 




Avg. 1460 

138.5 

7D18 

II 

10^ 

1682 

125 

7D19 

II 

10^ 

1650 

130 

7D20 

It 

10^ 

1378 

129 


Avg . 1570 

12S 


FF Fatigue Failure 

*Residual Strength Tests Conducted After N Cycles 
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Table 19 . 


Uniaxial Compression/Compression Lifetime 
and Residual Strength/Stiffness Data for 
[0/90/0] 2 g T-300/5208 Laminate (f=30Hz , R=10) 



Residual 

Strength 

Residual 
Modulus 
E (GPa) 

X 

Poisson ' s 
Ratio 

V 

xy 

Specimen 

S 

N, Cycles 

(MPa) 

X 

5BC5 

0 . 8 

150 




SBC 4 

II 

640 




5AL1 

II 

15,800 




SAL 4 

II 

10 




5BC6 

It 

2,770 




5AL2 

II 

1,110 




5BL4 

0.67 

5,920 




SBC 2 

II 

25,290 


- 


5AL6 

II 

6* 

5x10 

816 

97 

0.172 

5CL1 

II 

2,009,050 

1061 

79 

0.05 

5BL2 

II 

5x10^ 

966 

93 

0.03 

SAL 3 

II 

2,470 

- 

- 

- 




Avg . 1014 

86 

0.04 

5BL6 

0.5 

6* 

5x10 

829 

83 

0.183 

SAC 3 

II 

6* 

5x10 

872 

87 

0.06 

5BL3 

II 

5x10^* 

988 

§1 

0.092 




Avg . 897 

86 

0.112 

5AL1 

II 

10^ 

852 

97 

0.072 

5BL5 

II 

10^ 

865 

9 3 

■ 0.044 

SAC 2 

II 

10® 

842 

90_ 

0.109 




853 

93 

0.075 


*Runout 











Table Bl. Design of a Partially Notched [±45./0.]^ 
Laminate to Assure Failure due to J 

Interlaminar Stress t 

xz 


% of 0° Layers 

40% 

60% 

80% 

Initial Failure Stress 
for Laminate (MPa) 

(±45° Layers in Shear) 

393.0 

554 . 7 

730 . 8 

Ultimate Failure Stress 
for Laminate from 
Netting Analysis (MPa) 
(0° Layer in Tension) 

579.2 

868 . 7 

1158.3 

Failure Stress of ±45° 
Laminate in the Notched 
Area (MPa) 

124.1 

124 . 1 

124 . 1 

Corresponding Force 

F±^5.(N) 

(Specimen Width=25 . 4mm) 

2642.2 

1761.5 

880.7 

Value of the Force , F , 

s 

when Failure Occurs by 
Interlaminar Shear (N) 

( for max. T =68.9 5MPa) 

xz 

2139 . 6 

1347.8 

653 . 9 

F /F 

±45°'^ s 

1.23 

1 . 31* 

1.35 


*The 60% configuration is chosen for the experiments. 
Note: All failure stresses are for uniaxial tension. 

Lamina Thickness = 0.14 m.m 
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Figure 1. Minimechanics Model for a Notched Laminate 







A transverse 
crack through 
the thickness 
(catastrophic ) 




An axial crack 
in any layer 
between the core 
and overstressed 
regions 

(could grow in a 
stable manner till 
failure) 



o 

X 


A delaminat ion between 
any two layers in the 
core region 







Figure 4 . Damage Zone Sequence and Local Coordinates 






ELASTIC ANALYSIS 
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Figure 5. Various Damage Growth Combinations in a Notched Laminate 























virgin Data Input 


Elastic Analysis 



Diagnose inplane 
and/or interlaminar 
plastic damage 


Continue analysis with damage 


'ransverse Laminate Failure ? 


Diagnose a second inelastic 
damage zone, or a debond 
and/or inplane crack 


Continue analysis till 
crack/debond grows in an 
unstable fashion (axial 
laminate failure) , or till 
transverse laminate failure 


Stop 


Input properties after N cycles 


Figure 6 . Flow Chart for Failure Mode 
and Strength Predictions. 
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Laminate 

Displacement 


Figure 7 . Transverse Displacement Distributions .Near Sp^ 
Edge, [45/0/"45]^ Laminate (ref. 4 ), Showing 

Average Lamina Displacements, u^ , First Inter 
Shear Strain, y , and Determination of Parame 



Figure 8 . 


Elastic-Perf ectly Plastic Approxinaation 



Different Linear Approximations of the In-Plane and 
Interlaminar Shear Strain Curve 
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- — ■ — ► Inplane plasticity in the ±45° 

layers ( 0 - 2 ) 

Inplane plasticity in the 0° layers (a^^) 


[0„/±45/0] Laminate 
I s 


Trial Data Input 






Inpleine crack in 
the 0° layers (r. ) 


Inplane plasticity in the 
±45° layers (a^) 


Iriplane plasticity in the 
0° and 90° layers 


Trial Data Input 


[0/±45/90]__ Laminate 
s 


a (cm) 


Figure 10. Variation of Applied Laminate Force with 
the Total Axial Damage Zone Size 





Transverse Failure 


[a = 0.635 cm) 


Plasticity in the 
0° plies 


Transverse 

Failure 


(a = 1.27 cm) 


Plasticity in the 0° plies 


(a = 1.91 cm) 


Transverse 

Failure 


“Plasticity in the 0° plies 


[±45/0/0]^ Laminate 


Actual Data Input 


□ Transverse Failure 

(Secant Modulus Approach) 


a (cm) 


Figure 12 


Variation of Applied Laminate Force with 
the Total Axial Damage Zone Size 






























y = 0.281 


y = 0.625 


[0/ + 45/--45J Laminate 


Uniform value of a at y = 0.281 

yz 

from the approximate analysis 


( I Uniform value of at y = 0.625 

\y^ from the approximate analysis 


-•—Edge 
of the notch 




0.281 


0.625 


Figure 19 . Average for Unit Applied Average Stress 

at the Interface Between 0° and 45° Layers 
for Different Values of y. 









0 .08 


Circular Hole 
X = 0.125 (Edge of Hole) 
Order of Integration = 2 


0.06 

a 

yz 

0.04 

0.02 

0.00 



0.0 0.4 0.8 1.2 1.6 


Y 


Variation of o with y at the 0/+45 
y z 

Interface in a Notched [0/+45/-45] Laminate 


Figure 22 . 




0.07 


Circular Hole 
X = 0.125 (Edge of Hole) 
Order of Integration = 4 



• \J s. 

0.0 0,4 0.8 1.2 1 

y 


Figure 23 . Variation of with y at the 0/+45 

Interface in a Notched [0/+45/-45] 
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Figure 29. Plastic Elements in a [0/90] g Notched Laminate for P=28.06E3 
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Figure 31. 




Cel 1 



Figure 32. 

Loading Arrangement 
for the Modified 
IITRI Compression 
Test Fixture 




Actuator 


10 3 













Figure 35. [±45]^ T-300/5208 Laminate Failure Modes for Uniaxial 

Tension for Two Different Specimen Geometries 



Fiqure 36. [±45,^] T-300/5208 Laminate Compression Failures Showing 

lb S 

Microbuckling of Laminae 






















Ill 


[+45/0-] T-300/5208 Laminate 












Interlaminar Debond 



Figure 42 . 


Photomicrograph of a Failed [±45/02] ^ 
Partially Notched Laminate (xlOO) 
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90 



I Dnpissy 


Cycles 

Figure 45 . Residual Strength - S Data for tOl T-300/5208 Laminate 






T300/5208 

[o/eo/olgg 

Compression Fatigue 



Figure 47 . S-N Data for [0/90/0] 2 g T-300/5208 Laminate 





Figure 48. Residual Strength - N Data for [0/90/0] 2 ^ T-300/5208 Laminate 
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Figure Partially Notched [±45i/0j] Laminate 

for Obtaining Interlaminar ®Shear 
Stress (t ) Data (not to scale) 
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